Color confinement and dual superconductivity in full QCD 
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I. INTRODUCTION 

In a series of papers which we shall refer to as 

I, II, III respectively, we have demonstrated by numerical 
simulations that quenched QCD vacuum is a dual su- 
perconductor in the confining phase, and goes to normal 
state at the deconfinement transition (see also Ref. [Q). 
More precisely we have constructed an operator fj, car- 
rying magnetic charge and we have measured its vac- 
uum expectation value (/i). For T < Tc (/i) =/= 0, for 
T >Tc in) = and, approaching T^, (fi) ~ (1 - T/T^Y 
{6 — 0.50(3)). Magnetic charge is defined by a proce- 
dure called Abelian projection: it associates {Nc — 1) 
U{1) magnetic symmetries to any operator in the ad- 
joint representation j5|. A priori magnetic symmetries 
corresponding to different Abelian projections (different 
choices of cj)) are independent. In I, II, III we have shown 
that the behaviour of (/i) , including the value of the index 
(5, is independent of the Abelian projection. 

There is general agreement on the order-disorder na- 
ture of the deconfining transition in the quenched case. 
The popular order parameter is the Polyakov line (L); 
the symmetry involved is Z^. Alternatively the dual ('t 
Hooft) line (L) can be used as a disorder parameter 
(order parameter of the disordered phase) corresponding 
to the dual symmetry. Our (/x) is also a good disor- 
der parameter, and in fact it coincides numerically with 

In full QCD, i.e. in the presence of dynamical quarks, 
the situation is less clear. Zn and Zn symmetries are 
explicitely broken by the very presence of the quarks. At 
zero quark mass there is a phase transition at some Tc 
involving chiral symmetry: for T < Tc chiral symmetry is 
spontaneously broken, the pseudoscalar octet being the 
Goldstone particles, and for T > Tc it is restored. Quark 
masses do break chiral symmetry explicitely. It is not 
clear theoretically what the chiral transition has to do 
with the deconfinement transition. However, the suscep- 



tibilities of different quantities (the Polyakov line (L) , the 
chiral condensate) have been measured, and all of them 
have a maximum at the same value Tdniq) for any value 
of rriq. Above a certain value of (m^ > 3 GeV) the 
transition is first order, as in the quenched case, and (L) 
still works as an order parameter. At niq ~ the transi- 
tion is presumably second order. At intermediate values 
the susceptibilities which have been considered show a 
maximum at Tc, but it does not go large at increasing 
volume. The indication is then that there is no transi- 
tion but only a crossover. 

A natural question is then if dual superconductivity 
is a symmetry for the transition in full QCD as it is in 
the quenched case. In the spirit of the Nc — *■ oo limit, 
one would expect that the mechanism of confinement be 
the same as in quenched QCD, the idea being that the 
structure of the theory is the same as that in the limit 
iVc ^ oo at g'^Nc = A fixed: at finite Nc small differences 
are expected with respect to the limiting case. Quark 
loops are nonleading in the expansion. The mechanism 
of confinement should be approximately Nc independent 
and the same with and without dynamical quarks. 

The disorder parameter (/x) can be constructed in full 
QCD exactly in the same way as in the quenched case 
(see Section 2). At a given temperature T, (/Lt) has to be 
computed in the infinite volume limit. We have investi- 
gated the region T < Tc, where we find 

hm (m) ^ (1) 

and T > Tc^ where we find 

hm (m) = (2) 

V^oo 

as will be shown in detail below. Notice that the limit Eq. 

is not within errors but exact. Indeed we measure, 
instead of (/x), the quantity p = ^ln(/L(), and we find 
that it tends to — oo as p = —kNs -\- k' {k > 0) as the 
spatial size of the sample N^ oo. 
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The finite size scaling analysis in the critical region is 
under investigation, to study the nature and the order of 
the transition. 



II. DISORDER PARAMETER 

The operator /i is defined in full QCD exactly in the 
same way as in the quenched theory 



Z = J {VU)e 
Z = I {VU) 



■PS 



(3) 



Z is obtained from Z by changing the action in the time 
slice xq, S ^ S = S + AS. In the Abelian projected 
gauge the plaquettes 

= Ui{x,xo)Uo{x + i,xq)U^(x,xo + 6)U^{x,xo) (4) 

are changed by substituting 

U,{x, xo) ^ C/,(x, xo) = U,{x, xo)e'^''^^'^'y^ (5) 

where b{x — y) is the vector potential of a monopolc con- 
figuration centered at y in the gauge V5 = 0, and T is 
the diagonal gauge group generator corresponding to the 
monopole species chosen. In SU(2) T = (73/2, in SU(3) 
T = A3/2 or {V3Xs - A3)/2. In the generic SV{N) case 
the procedure is explained in Ref. Unlike the Z^^ 
centre symmetry, the U{1) magnetic symmetry defined 
after Abelian projection is a good symmetry also in pres- 
ence of dynamical fermions. It can be shown that, as 
in the quenched case, /Li adds to any configuration the 
monopole configuration h{x — y). If the magnetic sym- 
metry is realized a la Wigner, (/i) = if /i carries non 
zero net magnetic charge. Then (/z) 7^ means Higgs 
breaking of the U(l) symmetry. Therefore (/i) can be 
a correct disorder parameter for the transition to dual 
superconductivity also in full QCD. 

III. NUMERICAL RESULTS 

We have measured (/i) with two flavours of degenerate 
staggered fermions on x 4 lattices, with different val- 
ues of Ns {Ns — 12,16,32) and of the bare quark mass niq. 
In particular we have chosen, in the transition region, to 
vary the temperature, T = 1/ {Nta{(3,mq)), moving in 
the niq) plane while keeping a fixed value of m^/mp. 
To do this and to extract the physical scale we have used 
fits to the TOp and m-n masses pubhshed in jl^]. We 
present here results obtained at m.^ /mp ~ 0.505: in this 



case, at Nt = 4, the f3 corresponding to the transition is 
approximately Pc ~ 5.35 |ll| . Preliminary results have 
been already presented in | [12| . 

Instead of (^) we measure the quantity 



It follows from Eq. 



|) that 



(6) 



(7) 



the subscript meaning the action by which the average is 
performed. In terms of p 



(fj.) = exp 



p(/3')d/3' 



(8) 



A drop of (fi) at the phase transition corresponds to a 
strong negative peak of p. 

We have used the R version of the HMC algorithm for 
our simulations p^ . Some technical complications arise 
in the computation of the second term on the right hand 
side of Eq. (^. In the evaluation of {S)g, C*-periodic 
boundary conditions in time direction have to be used 
for the gauge fields and this requires C* boundary condi- 
tions in temporal direction also for fermionic variables (in 
addition to the usual antiperiodic ones), in order to en- 
sure gauge invariance of the fermionic determinant. This 
implies relevant changes in the formulation and imple- 
mentation of the HMC algorithm which are explained in 
detail in Ref. Hi- 

We have chosen the Polyakov line as the local adjoint 
operator which defines the Abelian projection. Actu- 
ally, calling L{x,Xq) the Polyakov line starting at point 
(x, a;o), the Abelian projection is defined by the operator 
L{x, xo)L*(x, Xo), which transforms in the adjoint repre- 
sentation when using C* boundary conditions. 

The use of a modified gauge action also implies changes 
in the molecular dynamics equations. One has to main- 
tain the modified hamiltonian containing S constant. A 
change in any temporal link indeed induces a change in 
L(x, xq) and hence in the Abelian projection defining the 
monopole field. Therefore the dependence of S on tem- 
poral links is non trivial and the equations of motion for 
the temporal momenta become more complicated. 

Fig. 1 shows p for a 32"^ x 4 lattice, and the chiral 
condensate as a function of (3. The negative peak of p is 
clearly at the same value of (3 where (ipip) drops to zero. 

Fig. 2 shows the plot of p for different spatial sizes Ng. 
For larger lattices the peak becomes deeper and the value 
of p at high P lower. 
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FIG. 1. Chiral condensate (open circles) and p (filled cir- 
cles) on the 32^ x 4 lattice. 



FIG. 2. Behaviour of p around the phase transition at var- 
ious lattice sizes. 



An analysis of p at large /3's as a function of is 
shown in Fig. 3, for different masses of the staggered 
fermions used in the simulations. For net magnetic 
charge ^ 



p ~ ~kN, + k' {k> 0) 



(9) 



and is practically independent of the quark mass within 
errors. For net charge zero (e.g. nionopole-antimonopole 
pair) p stays constant at large Ns . Going back to Eq. (||) 
this means that (p) is strictly zero in the infinite volume 
limit for non zero magnetic charge, and can be ^ for ex- 
citations with zero net magnetic charge. This statement 
is based on the analysis of many different excitations with 
different magnetic charges, and Fig. 3 is only an example. 
The magnetic symmetry is therefore realized d la Wigner 
for T > Tc and the Hilbert space is superselected. Notice 
that: 

(1) (p) can only be strictly zero in the infinite volume 
limit (Lee- Yang theorem ||l5||). 

(2) If we were measuring (p) directly we would find 
zero within large errors. Looking instead at p we can 
unambigously check Eq. (|^), which means that (p) is 
strictly zero as — s- oo. 

For T < Tc (p) 7^ if p stays constant and finite 
with increasing volume. This is what indeed happens 
as shown in Fig. 4. Nothing spectacular can happen at 
larger volumes, since no larger length scale exists in the 
system. 

Around Tc a finite size scaling analysis is required to 
get information on the order of the transition as well 
as to measure the critical indices. The problem is more 
complicated than in the quenched case, since an extra 
scale, the quark mass, is present. The program is on the 
way on a set of APEmille machines. Some qualitative 
features are shown in Fig. 2. 
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FIG. 3. Weak coupling behaviour of p at various lattice 
sizes. 



IV. CONCLUSIONS 

The preliminary data reported in this paper contain 
enough information to state that dual superconductiv- 
ity is at work as a confinement mechanism in QCD with 
dynamical quarks, in the same way as in the quenched 
theory [I,II,III]. For T > Tc the Hilbert space is supers- 
elected with respect to magnetic charge, for T < Tc the 
symmetry is Higgs broken. 

Dependence of the disorder parameter on the choice of 
the Abelian projection and nature of the transition are 
under investigation. 
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FIG. 4. Strong coupling behaviour of p at various lattice 
sizes and am = 0.1335. 
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